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Abstract.  Recent  work  has  developed  a  new  framework  for  the  impact  of  vertical  wind  shear  on  the 
intensity  evolution  of  tropical  cyclones.  A  focus  of  this  framework  is  on  the  frustration  of  the  tropical 
cyclone’s  power  machine  by  shear-induced,  persistent  downdrafts  that  flush  relatively  cool  and  dry 
(lower  equivalent  potential  temperature,  6e)  air  into  the  storm’s  inflow  layer.  These  previous  results 
5  have  been  based  on  idealised  numerical  experiments  for  which  we  have  deliberately  chosen  a  simple 
set  of  physical  parameterisations.  Before  efforts  are  undertaken  to  test  the  proposed  framework 
with  real  atmospheric  data,  we  here  survey  and  diagnose  five  additional  numerical  experiments  with 
some  modifications  of  the  experimental  setup  to  assess  the  robustness  of  our  previous  results.  The 
modifications  comprise  the  values  of  the  exchange  coefficients  of  surface  heat  and  momentum  fluxes, 
1 0  the  inclusion  of  experiments  with  ice  microphysics,  and  the  consideration  of  weaker,  but  still  mature 
tropical  cyclones. 

In  all  experiments,  the  depression  of  the  inflow  layer  9e  values  is  significant  and  all  tropical 
cyclones  exhibit  the  same  general  structural  changes  when  interacting  with  the  imposed  vertical  wind 
shear.  Tropical  cyclones  with  a  higher  downdraft  activity  exhibit  a  more  pronounced  depression  of 
15  inflow  layer  9e  outside  of  the  eyewall  in  our  experiments.  The  magnitude  of  the  9e  depression 
underneath  the  eyewall  early  after  shear  is  imposed  in  our  experiments  correlates  well  with  the 
magnitude  of  the  ensuing  weakening  of  the  respective  tropical  cyclone.  Based  on  the  evidence 
presented,  it  is  concluded  that  the  newly  proposed  framework  is  a  robust  description  of  intensity 
modification  in  our  suite  of  experiments. 
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1  Introduction 


1.1  A  new  framework  for  intensity  modification  of  tropical  cyclones  in  vertical  wind  shear 

The  intensity  evolution  of  tropical  cyclones  (TCs)  embedded  in  an  environmental  flow  with  vertical 
wind  shear  poses  an  important  forecast  problem.  To  date,  the  physical  processes  that  govern  this 
intensity  evolution  are  not  well  understood.  In  a  recent  publication  (Riemer  et  al.,  2010,  referred  to 
25  as  RMN10  hereafter)  a  new  framework  for  the  intensity  modification  of  TCs  in  vertical  wind  shear 
was  proposed.  While  previous  hypotheses  have  focused  almost  exclusively  on  processes  above  the 
boundary  layer,  RMN10  showed  that  vertical  wind  shear  has  a  significant  impact  on  the  thermody¬ 
namic  properties  of  the  TC’s  inflow  layer1.  Strong  and  persistent,  shear-induced  downdrafts  flush 
the  inflow  layer  with  low  moist-entropy  (or  equivalent  potential  temperature,  9e)  air  -  “anti-fuel”  for 
30  the  TC  power  machine.  The  replenishment  of  this  air  by  surface  heat  fluxes  is  not  complete  while  the 
air  spirals  towards  the  storm  centre  leading  to  a  reduction  of  9e  values  within  the  inner-core  updrafts 
in  the  azimuthal  mean.  When  viewed  from  the  simplified  perspective  of  an  idealised  Carnot-cycle  a 
decrease  in  the  potential  intensity  of  the  TC  can  be  expected  (Tang  and  Emanuel,  2010). 

The  downdrafts  that  flush  the  inflow  layer  with  low-0e  air  are  associated  with  a  quasi-stationary, 
35  azimuthal  wave-number  1,  convective  asymmetry  outside  of  the  eyewall.  This  convective  asymmetry 
is  reminiscent  of  the  ’stationary  band  complex’  (Willoughby  et  al.,  1984)  and  for  this  reason  will  be 
referred  to  as  the  SBC  hereafter.  The  SBC  and  the  associated  downdraft  pattern  extends  outwards 
to  approximately  150  km  -  180  km.  On  account  of  the  swirling  winds,  the  updrafts  within  the  SBC 
rise  along  a  helical  path.  A  simple  kinematic  model  (Riemer  and  Montgomery,  2011)  suggests 
40  that,  at  low  levels  above  the  inflow  layer,  environmental  low -9e  air  approaches  the  storm  closely 
underneath  the  helical  updrafts  in  the  downshear  and  downshear-left  quadrants.  Evaporative  cooling 
of  precipitation  falling  out  of  the  SBC,  and  into  this  unsaturated  air,  most  likely  drives  the  strong  and 
persistent,  vortex-scale  downdrafts  that  cool  and  dry  the  TC’s  inflow  layer. 

In  RMN10,  evidence  was  presented  that  the  formation  of  the  SBC  is  connected  to  the  tilt2  of  the 
45  TC  vortex.  When  subject  to  vertical  wind  shear,  the  tilt  evolution  of  the  TC  vortex  was  shown  to  be 
governed  to  zero  order  by  asymmetric-balance  dynamics.  Vortex  tilt  can  thus  be  described  in  terms  of 
a  standing  vortex-Rossby-wave  (VRW)  pattern  with  azimuthal  and  vertical  wave-number  1  structure. 
Even  though  the  tilt  of  the  inner  core  is  only  10km-20km,  the  VRW  pattern  associated  with  the  tilt 
of  the  outer  vortex  extends  outwards  to  radii  of  150  km-200  km  (Fig.  15a  in  RMN 10)3.  After  vertical 

'in  the  current  study,  the  term  'inflow  layer’  will  refer  to  the  layer  of  strong  inflow  associated  with  surface  friction  in  the 
lowest  1  km  -1.5  km  where  the  departure  from  gradient  wind  balance  in  the  radial  direction  is  generally  a  maximum  (e.g., 
Bui  et  al.  (2009),  their  Figs.  5  and  6  and  accompanying  discussion). 

2  In  RMN  10,  the  tilt  is  defined  as  the  vector  difference  between  the  location  of  the  vorticity  centroids  at  10  km  and  1  km 
height. 

3  A  considerable  tilt  of  the  outer  parts  of  TC-like  vortices  with  small  tilt  of  the  inner  core  has  been  noted  previously  by 
other  authors  also  (see  e.g.  Fig.  6a  in  Jones  (1995),  Fig.  10  in  Reasor  et  al.  (2004)). 


2 


50  wind  shear  is  suddenly  imposed,  the  TC  settles  quickly  (approx,  within  3  h)  into  a  quasi-equilibrium 
tilt  direction  that  is  to  the  downshear-left.  Thus,  the  positive  vorticity  anomaly  associated  with  the 
standing  wave-number  1  VRW  pattern  is  found  downshear-left  at  upper  levels  and  downshear-right 
at  low  levels.  Frictional  convergence  associated  with  this  positive  low-level  vorticity  anomaly  then 
provides  forcing  for  the  asymmetric  convection  within  the  SBC.  The  foregoing  results  suggest  an 
55  important  connection  between  the  asymmetric -balance  dynamics  governing  the  tilt  evolution  of  the 
vortex  subject  to  vertical  shear  and  the  thermodynamic  impact  on  the  inflow  layer. 

The  findings  of  RMN10  were  derived  from  idealised  numerical  experiments  following  the  setup 
of  Bender  (1997)  and  Frank  and  Ritchie  (2001).  In  this  experimental  setup,  a  uni-directional  vertical 
wind  shear  profile  is  imposed  suddenly  on  a  mature  model  TC.  In  an  attempt  to  minimise  complex- 
60  ity  and  to  understand  the  fundamental  processes,  RMN10  deliberately  chose  a  simple  set  of  model 
parameterisations.  As  is  the  nature  of  results  derived  from  such  an  idealised  approach,  uncertainties 
exist  whether  the  newly  proposed  pathway  operates  in  the  real  atmosphere.  Although  some  prelimi¬ 
nary  observational  support  was  offered  in  RMN10  a  more  comprehensive  test  of  the  new  hypothesis 
with  atmospheric  data  awaits  future  studies.  Before  undertaking  such  efforts  we  here  survey  and 
65  diagnose  five  additional  numerical  experiments  with  some  modifications  of  the  experimental  setup 
to  assess  the  robustness  of  RMNlO’s  results. 

1.2  Summary  of  numerical  experiments 

The  numerical  experiments  in  RMN10  employed  a  microphysics  scheme  that  considers  “warm 
cloud"  processes  associated  with  the  collision-coalescence  process  only  (hereafter  termed  “warm 
70  rain”).  The  current  study  considers  experiments  with  ice  microphysics  also.  Adding  ice  processes  in¬ 
troduces  the  latent  heat  of  fusion/  sublimation,  new  classes  of  hydrometeors  with  various  fall  speeds, 
and  additional  pathways  to  form  precipitation.  Details  of  the  employed  ice  microphysic  scheme  are 
given  in  section  2.2.  Arguably,  the  most  significant  of  the  above  processes  is  the  additional  conver¬ 
sion  of  latent  heat  associated  with  the  freezing  and  melting  of  ice.  Ice  microphysics  produce  stronger 
75  updrafts  than  a  warm-rain  scheme  due  to  the  additional  latent  heat  release  associated  with  freezing. 
In  addition,  the  inclusion  of  precipitation  in  the  form  of  ice  may  significantly  enhance  the  formation 
and  intensity  of  downdrafts  (Srivastava,  1987). 

Previous  studies  have  shown  that  the  inclusion  of  ice  microphysics  in  numerical  models  may 
alter  the  radial  profile  of  the  swirling  winds  of  TCs  as  compared  to  a  warm-rain  experiment  (e.g. 
80  Fovell  et  ah,  2009).  In  general,  it  is  not  well  understood  how  these  differences  arise.  Fovell  et  al. 
have  proposed  a  mechanism  that  is  based  on  the  sedimentation  of  ice  particles  from  the  outflow 
anvil,  the  subsequent  sublimentation  of  these  ice  particles,  and  radiative  feedback.  In  our  suite 
of  experiments,  however,  radiative  processes  are  neglected.  Nevertheless,  potentially  significant 
differences  in  the  TC-vortex  structure  between  the  warm-rain  and  the  ice  microphysic  experiments 
85  do  arise,  indicating  that  the  mechanism  suggested  by  Fovell  et  al.  is  not  of  primary  importance.  The 
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potentially  important  differences  in  the  TC-vortex  structure  and  their  implication  for  the  interaction 
with  vertical  wind  shear  will  be  discussed  in  some  detail  in  section  3.1. 

Recent  studies  have  highlighted  the  importance  of  the  representation  of  the  planetary  boundary 
layer  in  numerical  models  for  vortex  intensification  and  structure  changes  of  simulated  storms  (Braun 
90  and  Tao,  2000;  Smith  and  Thomsen,  2010;  Montgomery  et  al.,  2010;  Smith  et  al.,  2011).  The  ex¬ 
change  coefficients  for  enthalpy,  Ck,  and  momentum,  6'/),  play  a  most  prominent  role.  They  govern 
the  energy  extraction  from  the  ocean  surface  and  the  dissipation  of  kinetic  energy  in  the  frictional 
boundary  layer,  respectively.  In  the  framework  of  RMN10,  the  values  of  Ck  and  C d  can  be  expected 
to  play  an  important  role  in  the  replenishment  of  the  depressed  boundary  layer  0e  air  while  spiralling 
95  towards  the  inner-core  updrafts.  The  replenishment  rate  is  proportional  to  Ck4-  The  strength  of  the 
inflow  within  the  frictional  boundary  layer  is  determined  in  part  by  Cd-  Smaller  values  of  Co  may 
lead  to  less  radial  inflow  and  thus  to  an  increase  of  the  replenishment  time.  Assuming  that  Ck  is  not 
changed,  the  replenishment  of  depressed  inflow  layer  9e  values  is  then  more  complete  than  for  larger 
values  of  Cd.  Based  on  these  considerations,  one  would  expect  that  the  shear-induced  flushing  of 
1 00  the  inflow  layer  with  low -6e  air  has  a  more  pronounced  impact  on  TC  intensity  for  smaller  values  of 
Ck  and  larger  values  of  Cd-  The  interpretation  of  the  importance  of  Ck  and  C-d  in  an  axisymmetric 
theory  for  a  TC  in  vertical  wind  shear  by  Tang  and  Emanuel  (2010)  is  consistent  with  the  foregoing 
considerations. 

Considerable  uncertainties  exist  about  the  values  of  Ck  and  Cd  at  the  high  wind  speeds  occurring 
105  in  TCs.  Deacon’s  formula  for  C-d  was  employed  in  RMN10  and  the  ratio  of  Ck/Cd  was  set 
to  unity.  Recent  results  derived  from  observational  data  (Black  et  al.,  2007;  Zhang  et  al.,  2008) 
strongly  indicate  that  Deacon’s  formula  overestimates  Cd  at  high  wind  speeds  and  that  Ck/Cd 
is  considerably  less  than  unity.  The  observational  data  exhibit  considerable  scatter  and  are  limited 
to  wind  speeds  of  marginal  hurricane  strength5.  To  address  the  likely  overestimation  of  Ck  and 
1 10  Ck/C-d  in  RMN10,  we  employ  in  this  study  values  of  Ck  and  Ck/Cd  that  are  consistent  with, 
but  rather  on  the  lower  end  of  the  range  of  values  observed  by  Black  et  al.  and  Zhang  et  al..  Further 
details  are  provided  in  section  2.3. 

The  TCs  considered  in  RMN10  were  intense  and  rapidly  intensifying  at  the  time  when  the  vertical 
wind  shear  profile  was  imposed.  In  the  current  study  we  will  consider,  inter  alia,  weaker,  but  still 
115  mature  TCs.  Finally,  our  additional  experiments  now  exhibit  various  rates  of  intensification  just 
before  vertical  wind  shear  is  imposed.  We  did  not,  however,  find  a  sensitivity  of  the  subsequent  TC 
evolution  in  vertical  shear  on  this  preceding  intensification  rate  in  our  experiments.  Thus,  the  various 
rates  of  intensification  just  before  vertical  wind  shear  is  imposed  will  not  be  considered  hereafter. 

The  next  section  provides  an  overview  of  the  numerical  model  employed  and  the  additional  ex- 
1 20  periments  performed,  an  expanded  description  of  the  ice  microphysics  scheme,  and  the  values  for 

4The  replenishment  rate  is  proportional  to  the  air-sea  thermodynamic  disequilibrium  also. 

5  Some  significant  strides  towards  determining  these  coefficients  at  major  hurricane  wind  conditions  have  been  made 
recently  by  Bell  (2010)  as  part  of  his  Ph.D.  dissertation. 
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the  surface  exchange  coefficients.  In  section  3  we  compare  the  evolution  of  the  quintessential  in¬ 
gredients  of  RMNlO’s  newly  proposed  framework  in  the  different  numerical  experiments.  One 
experiment  with  ice  microphysics  is  examined  in  further  detail  in  section  4,  including  an  estimate  of 
the  spindown  time-scale  of  the  simulated  TC-vortex  based  on  simple  axisymmetric  theory.  Section 
1 25  5  contains  our  conclusions. 

2  Model  setup  and  additional  experiments 

2.1  A  brief  model  description 

For  the  numerical  experiments  we  employ  the  Regional  Atmospheric  Modeling  System  (RAMS), 
developed  at  Colorado  State  University  (Pielke  and  Coauthors,  1992;  Cotton  and  Coauthors,  2003). 
1 30  The  RAMS  is  a  three-dimensional,  non-hydrostatic  numerical  modeling  system  comprising  time- 
dependent  equations  for  velocity,  non-dimensional  pressure  perturbation,  ice-liquid  water  potential 
temperature  (Tripoli  and  Cotton,  1981),  and  mixing  ratios  of  total  water  and  hydrometeors. 

A  standard  first-order  turbulence  scheme  developed  by  Smagorinsky  (1963)  is  used  for  sub-grid- 
scale  mixing,  with  modifications  by  Lilly  (1962)  and  Hill  (1974)  that  enhance  diffusion  in  unstable 
135  conditions  and  reduce  diffusion  in  stable  conditions.  For  simplicity,  radiative  processes  are  ne¬ 
glected.  The  bulk  aerodynamic  formulas  are  used  to  calculate  momentum  and  sensible  and  latent 
heat  fluxes  from  the  sea  surface.  The  values  of  the  exchange  coefficients  are  given  in  subsection  2.3. 
We  either  employ  a  warm-rain  microphysics  scheme  based  on  the  scheme  of  Kessler  (1969)  or  the 
RAMS  standard  ice  microphysics  scheme  (subsection  2.2).  A  comprehensive  description  of  RAMS, 
1 40  the  surface  interaction  parameterisation,  and  the  warm-rain  scheme  is  given  in  RMN10. 

2.2  Ice  microphysics  in  RAMS 

The  RAMS  ice  microphysics  scheme  has  eight  water  categories;  vapour,  cloud  droplets,  rain,  pristine 
ice,  snow,  aggregates,  graupel  and  hail  (Walko  et  al.,  1995;  Meyers  et  ah,  1997).  The  liquid  phase 
categories,  cloud  droplets  and  rain,  may  be  supercooled.  Pristine  ice,  snow  and  aggregates  are 
145  assumed  to  be  completely  frozen,  whereas  graupel  and  hail  are  mixed  phase  categories,  capable 
of  consisting  of  ice  only,  or  a  mixture  of  ice  and  liquid.  All  hydrometeor  categories  except  cloud 
droplets  are  assumed  large  enough  to  fall.  Cloud  droplets  and  pristine  ice  are  the  only  categories 
to  nucleate  from  vapour.  All  other  categories  form  from  existing  hydrometeors,  but  once  formed 
may  also  grow  by  vapour  deposition.  Pristine  ice  crystals  form  from  heterogeneous  nucleation, 
1 50  homogeneous  nucleation,  and  secondary  production  that  involves  splintering  of  riming  crystals.  The 
snow  category  is  defined  in  RAMS  as  relatively  large  ice  crystals,  which  have  grown  from  the  pristine 
ice  category  due  to  vapour  deposition  and  may  continue  to  grow  by  riming.  Aggregates  are  defined 
as  ice  particles  that  have  formed  by  collision  and  coalescence  of  pristine  ice,  snow  and/or  aggregates. 
Snow  and  aggregates  retain  their  identity  with  moderate  amounts  of  riming,  but  convert  to  graupel 


5 


155  if  the  amount  of  riming  is  large.  In  addition  to  forming  by  heavy  riming  graupel  can  also  form  by 
partial  melting  of  snow  or  aggregates.  Graupel  is  allowed  to  carry  up  to  half  its  mass  in  liquid.  If 
the  percentage  becomes  larger,  by  either  melting,  riming  or  collection  by  rain,  a  graupel  particle  is 
re-categorised  as  hail.  Hail  is  formed  by  freezing  of  rain  drops  or  by  riming  or  partial  melting  of 
graupel.  Hail  is  allowed  to  be  a  large  fraction  liquid  and  converts  to  rain  when  completely  liquid. 

1 60  The  hydrometeor  size  distribution  in  RAMS  is  represented  by  a  generalised  Gamma  function.  A 
user  specified  shape  parameter  controls  the  width  of  the  distribution.  Also,  a  scaling  diameter  is 
often  specified  for  a  one-moment  scheme  that  is  related  to  the  mean  diameter  of  the  distribution, 
with  the  particular  relation  depending  on  the  chosen  spectral  parameters.  For  this  study  a  one- 
moment  scheme  is  employed  for  all  categories  except  cloud  droplets  and  pristine  ice.  A  scaling 
1 65  diameter  is  specified  and  the  number  concentration  is  diagnosed  using  the  predicted  mixing  ratio  of 
the  hydrometeor  scheme. 

A  heat  budget  equation  is  formulated  for  rain,  graupel,  and  hail  hydrometeors.  This  allows  the 
temperature  of  a  hydrometeor  to  differ  from  that  of  the  surrounding  ambient  air  temperature  due  to 
latent  heat  release  or  absorption  in  the  hydrometeor,  and  sensible  heating  by  collisions  with  other 
170  hydrometeors.  The  temperature  of  a  hydrometeor  often  differs  substantially  from  that  of  the  air, 
which  can  significantly  influence  the  rates  of  heat  and  vapour  diffusion,  as  well  as  the  amount  of 
sensible  heat  transfer  that  occurs  in  coalescence  of  hydrometeors.  Collision  and  coalescence  of 
hydrometeors  is  governed  by  the  stochastic  equation  derived  byVerlinde  et  al.  (1990). 

Gravitational  settling  of  hydrometeors  causes  them  to  fall  relative  to  air.  Terminal  velocities  in- 
1 75  crease  as  a  power  of  the  droplet  diameter,  and  are  based  on  empirical  data.  Pristine  ice  has  very  small 
fall  speeds.  Snow  crystals  fall  faster  than  pristine  ice,  but  still  relatively  slowly.  Aggregates  are  large 
crystals  that  have  moderate  fall  speeds.  Graupel  and  hail  have  high  fall  speeds.  Sedimentation  in 
RAMS  uses  a  Lagrangian  scheme  to  transport  the  mixing  ratio  from  any  given  grid  cell  to  a  lower 
height  in  the  vertical  column,  so  that  a  proportion  is  transferred  to  the  lower  adjacent  grid  cell. 

1 80  2.3  Modification  of  surface  exchange  coefficients 

Deacon’s  formula  for  the  dependence  of  the  drag  exchange  coefficient  C n  on  wind  speed  is  given 
by 

CD  =  1.1  x  10-3  +  4  x  10— 5 1  V^| 1 ,  (1) 

with  | HI1  being  the  wind  speed  in  ms-1  at  the  lowest  model  level  (here:  49  m).  In  RMN10  we  used 
1 85  Deacon’s  formula  and  Ck  =  Cd ■  For  reference,  this  formula  furnishes  values  of  the  exchange  coef¬ 
ficients  of  1.5  x  10-3,  2.1  x  1CT3,  and  3.1  xlO-3  for  10  ms-1,  25  ms-1,  and  50  ms-1,  respectively. 

Zhang  et  al.  (2008)  found  that  ’The  exchange  coefficient  of  enthalpy  flux  shows  no  significant 
dependence  on  wind  speed  up  to  hurricane  force  with  a  value  of  1.16  x  10-3.’  and  that  ’The  average 


6 


ratio  of  the  CkICjj  values  is  0.63, We  therefore  use  these  values  and  set 
190  Ck  —  1.16  x  10-3  =  constant  (2) 

and 

CK/CD  =  0.63.  (3) 

Experiments  that  employ  these  values  are  referred  to  as  CBLAST  below.  Note  that  Deacon’s  formula 
and  Ck/Cd  =  1  furnishes  values  of  Ck  and  Co  that,  for  ll^1  =  50  ms-1,  are  approx.  2.7  and  4.2 
195  times  higher,  respectively,  than  the  corresponding  CBLAST  values. 

2.4  Summary  of  experiments 

A  summary  of  the  experiments  examined  in  this  study  is  given  in  table  1 .  We  refer  to  experiments  that 
employ  warm  rain  and  Deacon’s  formula  for  the  exchange  coefficients  (with  Ck /Cd  =  1)  as  RMN. 
The  CBLAST  experiments  employ  warm  rain  also.  Experiments  that  include  ice  microphysics  are 
200  referred  to  as  ICE.  ICE  experiments  use  Deacon’s  formula  with  Ck  / Cjj  =  1. 

The  shear  profile  employed  in  this  study  has  a  cosine  structure  in  the  vertical  with  zero  winds  at 
the  surface  and  maximum  easterly  winds  at  12  km  and  above.  We  will  follow  the  convention  of  the 
TC  community  to  define  vertical  wind  shear  as  the  vector  difference  between  the  winds  at  upper 
and  lower  levels.  Here,  we  use  the  vector  difference  between  12  km  and  the  surface.  The  shear 
205  magnitude  is  thus  given  in  ms-1.  Unless  otherwise  noted,  we  will  consider  TCs  in  vertical  wind 
shear  with  a  magnitude  of  15  ms-1.  The  shear  profile  is  imposed  at  a  high  intensity6  of  68  ms-1, 
after  the  TCs  have  spun  up  in  quiescent  environment.  These  experiments  are  referred  to  as  RMNgg, 
CBLAST68,  and  ICE68,  respectively.  The  intensity  at  the  time  when  shear  is  imposed  will  be  referred 
to  as  “initial  intensity”  below.  Experiment  RMNf;*  is  the  15mps  case  of  RMN10  and  constitutes  our 
210  reference  experiment. 

We  have  performed  experiments  with  a  weaker  initial  intensity  of  54  ms-1  also,  referred  to  as 
RMN54  andCBLAST54.  Furthermore,  an  ICE  experiment  has  been  performed  with  10ms_1  vertical 
shear  fioICEgs)  which  will  be  compared  briefly  to  the  lOmps  case  of  RMN10  (here:  loRMNgs). 

3  Assessing  the  robustness  of  the  new  framework 

215  We  shall  focus  on  the  key  features  of  RMN  10’s  new  framework  for  intensity  modification  in  vertical 
wind  shear:  the  vortex  tilt,  the  occurrence  and  location  of  the  SBC,  the  flushing  of  the  inflow  layer 
with  low-0e  air,  the  ensuing  depression  of  inflow  layer  9e,  and  the  associated  weakening  of  the  TC. 
The  general  characteristics  of  these  key  feature  will  be  found  to  be  consistent  in  all  experiments. 

6The  storm  intensity  metric  in  this  study,  as  well  as  in  RMN  10.  is  the  maximum  azimuthally  averaged  tangential  wind 
speed  at  1  km  height. 
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As  can  be  expected,  differences  arise  in  the  specific  characteristics  of  these  key  features.  For  some 
220  of  these  differences  it  is  very  difficult  to  establish  the  causal  link  to  the  changes  in  the  experimental 
setup.  Most  prominently,  differences  between  the  evolution  of  the  TCs  in  RMN(iX  and  rCE(i8  are 
partly  due  to  the  enhancement  of  downdrafts  by  ice  microphysics  (Sec.  3.7)  but  very  likely  also 
due  to  changes  in  the  radial  profile  of  the  swirling  winds  at  the  time  when  vertical  wind  shear  is 
imposed.  Disentangling  the  individual  contributions  requires  a  more  complete  understanding  of  the 
225  processes  at  work,  and  extended  numerical  experimentation  and  diagnostic  analyses,  preferably  in  an 
ensemble-based  framework.  A  more  comprehensive  examination  of  such  differences  is  thus  deferred 
to  future  research. 

3.1  Importance  of  the  radial  profile  of  the  swirling  winds 

3.1.1  Impact  on  the  resiliency  of  TC-like  vortices 

230  It  is  generally  believed  that  the  presence  of  vertical  wind  shear  increases  the  interaction  of  a  TC  with 
dry  environmental  air  (Riemer  and  Montgomery,  2011,  and  references  therein).  Based  on  a  simple 
kinematic  model,  Riemer  and  Montgomery  (2011)  argue  that  a  TC  with  a  broad  radial  profile  of  the 
swirling  winds  should  have  a  stronger  ability  to  isolate  itself  from  the  detrimental  interaction  with 
dry  environmental  air  than  a  smaller-scale  TC.  Furthermore,  the  radial  profile  of  the  swirling  winds 
235  has  a  profound  impact  on  the  dynamic  resiliency  of  TC-like  vortices,  i.e.  the  ability  to  maintain 
vertical  coherence  in  the  presence  of  differential  advection  by  the  vertical  shear  flow  (e.g.  Jones, 
1995;  Reasor  et  ah,  2004).  These  studies  investigated  the  resiliency  of  initially  barotropic,  dry 
vortices  and  are  reviewed  briefly  below. 

Jones  (1995)  examined  vortex  evolution  in  terms  of  the  formation  and  subsequent  interaction  of 
240  potential  vorticity  (PV)  anomalies.  PV  anomalies  were  defined  as  deviations  from  the  azimuthal 
mean.  Due  to  the  differential  advection  by  the  vertical  shear  flow,  PV  anomalies  representing  the 
tilt  of  the  vortex  form  at  lower  and  upper  levels.  These  anomalies  mutually  interact  through  the 
vertical  penetration  of  their  “induced”  quasi-horizontal  circulations,  causing  the  vortex  to  precess. 
This  precession  helps  the  vortex  to  counteract  the  differential  advection  by  the  vertical  wind  shear. 
245  In  general,  the  vertical  interaction  of  PV  anomalies  depends  on  the  Rossby  penetration  depth.  The 
penetration  depth  increases  with  the  strength  of  the  background  rotation  and  with  the  horizontal  scale 
of  the  interacting  anomalies.  Jones  noted  that  it  is  not  entirely  clear  how  to  quantify  the  penetration 
depth  for  a  tilted  TC-like  vortex  but  concluded  that  intense  and  broad  vortices  are  more  resilient  than 
their  weaker,  smaller-scale  counterparts. 

250  Reasor  et  al.  (2004)  re-visited  the  problem  of  vortex  resiliency  and  re-alignment  using  VRWs  con¬ 
cepts  and  models.  Their  results  are  in  qualitative  agreements  with  Jones’s  finding  that  vortex  (and 
environment)  parameters  that  increase  the  penetration  depth  improve  the  resistence  of  a  vortex  to 
tilting  by  an  imposed  vertical  shear  flow.  In  addition,  Reasor  et  al.  demonstrated  that  vortex  preces- 
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sion  is  not  the  only  element  in  the  resiliency  of  TC-like  vortices.  An  inviscid  damping  mechanism 
255  of  the  vortex-tilt  mode  exists,  associated  with  PV  mixing  at  the  critical  radius  of  the  near-discrete 
VRW  that  represents  the  vortex  tilt.  The  critical  radius  is  typically  located  in  the  “skirt”  region  of 
the  vortex,  i.e.  in  the  region  of  small  PV  values  outside  of  the  high-PV  values  in  the  inner  core  (ap¬ 
prox.  outside  50  km  -  70  km  in  Fig.  1(a)).  The  inviscid  damping  mechanism  depends  on  the  radial 
gradient  of  the  skirt  PV  at  the  critical  radius  (Schecter  and  Montgomery,  2003).  Inviscid  damping, 
260  and  thus  re-alignment  of  the  vortex,  occurs  when  the  radial  PV  gradient  in  this  region  is  negative. 
The  damping  rate  increases  with  the  magnitude  of  the  gradient.  For  a  positive  skirt-PV  gradient  a  tilt 
instability  may  occur.  Reasor  et  al.  ’s  results  clearly  demonstrate  that  vortex  resiliency  may  crucially 
depend  on  the  radial  profile  of  the  swirling  winds. 

The  radial  profile  of  the  swirling  winds  has  been  confirmed  to  play  a  role  in  the  intensity  evolution 
265  of  vertically  sheared  TCs  in  idealised,  moist  numerical  experiments  (Wong  and  Chan,  2004).  In 
these  experiments  the  weakening  of  smaller-scale  TCs  is  found  to  be  more  pronounced  than  that 
of  their  broader-scale  counterparts.  Wong  and  Chan  based  their  interpretation  of  this  difference 
in  the  intensity  evolution  on  the  results  of  Jones  (1995)  derived  from  dry  TC-like  vortices.  To  us, 
however,  it  is  not  clear  what  the  underlying  processes  are  that  govern  the  intensity  evolution  in  Wong 
270  and  Chan’s  experiments.  In  particular,  the  role  of  moist  processes  and/or  the  associated  secondary 
circulation  remain  unclear. 

In  the  real  atmosphere,  an  indication  of  a  higher  susceptibility  of  smaller-scale  TCs  to  vertical 
shear  has  been  found  by  DeMaria  (1996)  based  on  a  regression  analysis. 

3.1.2  Diagnosed  radial  vortex  structure  in  our  experiments 

275  Because  the  radial  profile  of  the  swirling  winds  may  play  a  prominent  role  in  the  evolution  of  TCs  in 
vertical  wind  shear,  we  begin  the  examination  of  our  experiments  with  the  comparison  of  the  vortex 
structure  at  the  time  when  vertical  wind  shear  is  imposed.  The  radial,  azimuthal  mean  profiles  of 
the  swirling  winds  at  this  time  are  illustrated  in  Fig.  1.  In  CBLAST(i8,  the  radius  of  maximum 
winds  (RMW)  is  found  at  40  km  (Fig.  1(a)).  The  vorticity  maximum  is  35xl0~4  s-1  at  a  radius 
280  of  20  km.  This  vortex  structure  is  very  similar  as  in  RMNT(l8  (shown  in  Fig.  4  in  RMN10)7.  For 
the  weaker  TC  in  CBLAST54,  the  vorticity  structure  outside  a  radius  of  35  km  is  almost  identical  to 
CBLAST88.  Inside  this  radius,  vorticity  is  virtually  constant  with  a  value  of  approx.  20xl0-4  s_1. 
The  corresponding  RMW  is  found  at  48  km.  The  initial  radial  structure  of  the  TC  in  RMN54  is  very 
similar  (not  shown). 

285  The  vorticity  structure  of  the  TC  in  ICE68  exhibits  some  notable  differences  as  compared  to  the 
warm -rain  experiments  CBLAST()8  and  RMN(i8.  The  vorticity  maximum  is  found  inside  a  radius  of 
10  km  with  values  of  55  x  1 0“  1  s_1,  i.e.  approx.  60  %  higher  than  in  CBLAST(i8.  Radially  outside 

7  A  potential  impact  of  barotropic  instability  associated  with  the  reversal  of  the  radial  vorticity  gradient  in  the  inner  core 
on  the  intensity  evolution  of  the  sheared  TCs  in  our  experiments  is  arguably  negligible  (see  Sec.  3.1  in  RMN10). 
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of  this  maximum,  the  high  vorticity  values  decrease  rapidly.  Corresponding  to  the  smaller  radius  of 
maximum  vorticity,  the  RMW  is  smaller  (located  at  30  km)  than  in  CBLASTgg  and  RMNfj8  also. 
290  As  discussed  above,  the  resiliency  of  the  TC  vortex  may  crucially  depend  on  the  gradient  of  (po¬ 
tential)  vorticity  at  the  critical  radius  of  the  (near-discrete)  VRW  mode  that  represents  the  vortex  tilt. 
The  critical  radius  is  found  where  the  azimuthal  propagation  of  this  VRW,  i.e.  the  vortex  precession, 
equals  the  speed  of  the  swirling  winds.  The  precession  rate,  and  thus  the  critical  radius,  in  our  ex¬ 
periments  is  difficult  to  determine.  After  the  vortices  have  re-aligned,  a  short  precession  period  of 
295  2  h  -  3  h  is  indicated  in  all  experiments  (see  Fig.  4  below8).  Using  the  values  of  the  axisymmetric 

swirling  winds  at  1  km,  as  depicted  in  Fig.  1(a),  our  rough  estimate  of  the  critical  radius  is  then 
50  km  -  90  km. 

In  ICE,  the  vorticity  gradient  is  near-zero  between  50  km  and  65  km  while  in  CBLASTgg  and 
RMN.54  the  gradient  is  negative  and  of  large  magnitude  (Fig.  1(b)).  Between  65  km  and  100  km 
300  the  radial  vorticity  profiles  are  broadly  similar.  At  larger  radii  (not  shown)  the  average  vorticity 
gradient  is  positive  in  RMN54  (between  120  km  and  200  km)  and  ICE  (between  150  km  -  200  km). 
In  CBLAST,  in  contrast,  vorticity  decreases  monotonically  out  to  200  km  radius.  It  is  questionable, 
however,  whether  these  differences  at  a  radius  considerably  larger  than  our  estimate  for  the  critical 
radius  play  a  role  for  the  TC-vortex  evolution  in  our  experiments. 

305  From  the  discussion  in  this  subsection  it  is  clear  that  differences  in  the  evolution  of  the  TCs 
when  interacting  with  the  vertical  shear  flow  in  the  various  experiments  potentially  arise  due  to 
differences  in  the  radial  vorticity  profiles  before  shear  is  imposed.  The  weaker  vortices  in  RMN54 
and  CBLAST54,  as  well  as  the  smaller-scale  vortex  in  ICEfl8,  can  be  expected  to  be  more  susceptible 
to  tilting  than  their  stronger,  respectively  broader-scale  counterparts  in  RMNT(i8  and  CBLASTg*.  In 
310  ICEgg ,  the  susceptibility  to  tilting  is  potentially  enhanced  by  the  near-zero  gradient  of  the  skirt 
vorticity  inside  of  65  km.  In  Sec.  3.3  it  will  be  shown  that  the  tilt  evolution  is  broadly  consistent 
with  these  expectations.  Furthermore,  based  on  the  results  of  Wong  and  Chan  (2004),  one  can  expect 
a  more  pronounced  intensity  decrease  in  ICEfi8  based  on  the  radial  scale  of  the  vortex  alone.  The 
latter  point  needs  to  be  borne  in  mind  when  attempting  to  attribute  differences  between  ICE(j8  and 
315  RMNgg  to  their  different  representation  of  cloud  microphysical  processes. 

3.2  Intensity  evolution 

An  overview  of  the  intensity  evolution  in  our  experiments  is  presented  in  Fig.  2.  First,  a  brief  account 
of  the  intensity  evolution  in  a  quiescent  environment  is  given.  The  intensity  evolution  in  RMN  is 
described  in  some  detail  in  RMN10.  In  CBLAST  (Fig.  2(a)),  the  incipient  TC  slowly  intensifies  to 
320  25  ms”1  in  the  first  20  h.  A  period  of  rapid  intensification  ensues  during  which  the  TC  reaches  an 

intensity  of  53  ms”1  at  33  h.  In  the  subsequent  20  h  the  TC  consolidates  and  slightly  intensifies 

sNote  that  the  temporal  resolution  of  the  data  in  Fig.  4  is  1  h.  It  is  thus  most  likely  that  the  vortex  precession  is  undersam- 
pled. 
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to  56  ms-1.  Then,  the  TC  continues  to  intensify  at  a  moderate  rate  until  it  reaches  a  quasi-steady 
intensity  of  80  ms-1  -  85  ms-1  at  80  h.  In  ICE  (Fig.  2(b)),  the  TC  intensifies  to  approx.  40  ms-1 
at  30  h  after  a  short  gestation  period.  Very  rapid  intensification  ensues  and  the  intensity  reaches 
325  90  ms-1  at  50  h.  The  maximum,  quasi-steady  intensity  of  95  ms’  1  -  100  ms-1  is  reached  around 
72  h. 

In  all  experiments,  intensity  decreases  after  vertical  wind  shear  is  imposed.  As  in  RMN10,  all 
TCs  are  resilient  to  the  detrimental  impact  of  vertical  wind  shear.  After  a  period  of  weakening  the 
TCs  re-intensify;  particularly  rapidly  in  RMN549.  The  re-intensification  process  in  this  experimental 
330  configuration  is  discussed  in  some  detail  in  RMN10  (their  Sec.  6). 

To  facilitate  the  intercomparison  of  the  intensity  evolution  in  the  various  experiments,  the  intensity 
difference  to  the  respective  no-shear  experiment  is  presented  in  Figs.  2(c)  and  2(d).  Time  0  h  in 
these  figures  denotes  the  time  when  vertical  wind  shear  is  imposed  in  the  respective  experiment.  In 
all  warm-rain  experiments  (Fig.  2(c)),  the  characteristics  of  this  relative  weakening  are  similar  in 
335  the  first  20  h.  The  only  exception  being  that  the  weakening  in  CBLAST54  is  delayed  by  4  h.  This 
detail  will  be  examined  in  section  3.6.2.  After  20  h,  the  experiments  with  weaker  TCs  (RMN54  and 
CBLAST54)  continue  to  weaken.  In  contrast,  CBLAST(i8  maintains  its  intensity  difference  for  the 
next  20  h  and  RMN68  re-intensifies.  Based  on  these  plots  alone  it  might  therefore  be  concluded 
that  the  weaker  storms  are  more  susceptible  to  shear.  However,  this  is  not  necessarily  the  case  here. 
340  Note  that  20  h  after  shear  is  imposed  in  RMN54  and  CBLAST54  the  respective  reference  TCs  in 
quiescent  environment  are  still  rapidly  intensifying.  In  contrast,  the  reference  TCs  in  RMN(i8  and 
CBLASTes  intensify  only  slowly  20  h  after  vertical  wind  shear  is  imposed.  Thus,  the  differences  in 
the  relative  intensity  evolution  partly  arise  due  to  the  differences  in  the  contemporaneous  evolution 
of  the  reference  TC.  This  brief  discussion  illustrates  one  of  the  inevitable  challenges  to  interpreting 
345  different  TC  simulations  in  vertical  wind  shear. 

The  comparison  of  the  RMN  and  ICE  shear  experiments  is  less  complicated  because  the  intensity 
evolution  of  the  reference  runs  in  quiescent  environment  is  similar  (Figs.  2(b)  and  2(d)).  A  compar¬ 
ison  of  the  intensity  evolution  in  ioTCE(j8  and  ICE(l8,  respectively,  clearly  shows  that  the  TCs  in  the 
ICE  experiments  are  more  susceptible  to  shear  than  their  counterparts  in  the  warm-rain  experiments. 
350  The  differences  between  the  relative  weakening  in  xoICEgg  and  loRMNgg,  he.  the  experiments  with 
imposed  vertical  shear  of  10  ms-1,  is  on  average  approx.  10  ms-1  before  16  h,  and  on  average 
approx.  5ms”1  between  16  h  and  36  h  (Figs.  2(d)).  After  this  time,  it  is  no  longer  evident  that  the 
weakening  in  10ICE68  is  more  pronounced  than  in  10RMNg8. 

For  15  ms-1  shear,  the  difference  in  the  intensity  evolution  between  the  ice  and  warm-rain  ex- 
355  periments  is  striking.  The  maximum  relative  weakening  in  ICErl8  is  -60  ms-1  (Figs.  2(d)).  The 
TC  reaches  a  minimum  intensity  of  35  ms-1  (Figs.  2(b)).  For  RMN(i8,  the  respective  values  are  - 

^Because  the  focus  of  this  study  is  on  the  structural  changes  of  the  TCs  during  the  weakening  phase,  this  experiment  was 
discontinued  at  84  h. 
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25  ms-1  and  63  ms-1.  Notwithstanding  this  distinct  intensity  evolution,  it  will  be  shown  below  that 
the  stmctural  changes  in  ICE(j8  and  RMNgs  are  very  similar.  In  particular,  the  pronounced  weaken¬ 
ing  in  ICE68  is  associated  with  by  far  the  most  significant  depression  of  inflow  layer  (),.  (sections  3.6 
360  and  4.1).  The  evolution  in  ICEgg  is  thus  consistent  with  the  results  of  RMN10,  i.e.  the  shear-induced 
frustration  of  inflow  layer  9e  plays  a  crucial  role  in  the  intensity  evolution  in  ICE68.  Furthermore,  it 
will  be  argued  in  section  4.2  that  this  rapid  weakening  is  consistent  with  the  frictional  spindown  of 
an  axisymmetric  vortex  due  to  a  significant  reduction  of  the  inner-core  convective  mass  flux. 

3.3  Evolution  of  vortex  tilt 

365  Following  RMN10,  the  tilt  of  the  vortex  is  defined  here  as  the  vector  difference  between  the  vorticity 
centroids  at  10  km  and  1  km  height.  Details  on  the  calculation  of  the  vorticity  centroids  are  given  in 
RMN10. 

A  detailed  examination  of  the  tilt  evolution  in  RMi\’fj8  was  given  in  section  5.1  of  RMN10.  There, 
evidence  was  presented  that  the  tilt  evolution  is  governed  to  zero  order  by  asymmetric -balance 
370  dynamics.  The  TC  in  RMN(i8  settles  into  a  quasi-equilibrium  tilt  direction  to  the  downshear-left 
(south)  shortly  after  shear  is  imposed  (approx,  within  3  h).  A  quasi-equilibrium  tilt  direction  may  be 
achieved  when  the  cyclonic  precession  of  the  tilted  vortex  cancels  the  differential  advection  by  the 
vertical  shear  flow.  This  quasi-equilibrium  vortex  tilt  can  be  described  in  a  first  approximation  by  the 
balanced  dynamics  of  an  initially  barotropic,  dry  TC-like  vortex  (Jones,  1995;  Reasor  et  ah,  2004). 
375  In  RMN68,  the  tilt  amplitude  decreases  while  maintaining  the  quasi-equilibrium  tilt  direction.  The 
TC  virtually  re-aligns  within  approx.  24  h,  i.e.  after  this  time  the  tilt  amplitude  is  comparable  to  the 
tilt  amplitude  of  the  TC  in  the  no-shear  benchmark  experiment  (see  Fig.  5  in  RMN10).  Once  the  tilt 
amplitude  has  decreased  considerably,  a  high  variability  in  the  tilt  direction  ensues,  indicative  of  a 
high-frequency  precession  of  the  vortex. 

380  Figures  3  and  4  demonstrate  that  the  general  characteristics  of  the  tilt  evolution  are  similar  in  all 
of  our  experiments  in  the  first  24  h  -  30  h  after  vertical  wind  shear  is  imposed.  All  TCs  are  resilient 
during  this  time  and  the  tilt  direction  is  predominantly  to  the  South  (180°).  In  CBLAST(i8,  the 
maximum  tilt  magnitude,  as  well  as  the  temporal  tilt  evolution,  is  very  similar  to  that  of  RMN(i8. 
The  only  exception  being  that  the  local  maximum  at  23  h  is  more  pronounced  and,  thereafter,  the 
385  tilt  magnitude  remains  3  km  -  5  km  larger.  As  can  be  expected  from  the  discussion  in  Sec.  3.1,  the 
weaker  TCs  in  RMN54  and  CBFAST54  clearly  exhibit  a  larger  tilt  magnitude  in  the  first  24  h  than 
the  more  intense  TCs  in  RMN88,  CBLASTfj8,  and  rCE(j8.  Maximum  tilt  amplitudes  reach  40  km  in 
RMN54  and  30  km  in  CBFAST54.  Both  TCs  eventually  re-align. 

The  most  striking  feature  in  the  tilt  evolution  is  the  disruption  of  the  TC  vortex  after  30  h  in  ICE(l8 
390  (Fig.  3).  Until  this  time  the  TC  tends  to  re-align  slowly:  the  local  tilt  maximum  and  minimum 
at  25  h  and  30  h,  respectively,  are  smaller  than  their  respective  predecessors  at  9  h  and  15  h.  At 
30  h,  intensity  has  decreased  by  over  30  ms-1  as  compared  to  the  time  when  vertical  wind  shear  is 
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imposed  and  has  reached  the  smallest  absolute  value  in  all  experiments  (approx.  35  ms-1.  Fig.  2). 
Based  on  the  discussion  in  Sec.  3.1,  it  is  plausible  that  the  significantly  weakened  TC  can  no  longer 
395  withstand  the  differential  advection  associated  with  the  vertical  wind  shear.  Note  that  the  distinct 
weakening  in  the  first  20  h  precedes  the  disruption  of  the  TC  by  more  than  10  h.  After  30  h,  the 
TC  remains  as  a  coherent  entity  with  smaller  vertical  extent:  the  tilt  as  inferred  from  the  vorticity 
centroid  at  9  km  height  instead  of  10  km  height  does  not  exceed  35  km  (not  shown).  Arguably,  the 
survival  as  a  vertically  less  extensive  system  allows  for  the  rapid  re-intensification  of  the  TC  after 
400  the  inflow  layer  9e  values  have  recovered  (Sec.  4.1). 

The  tilt  magnitude  in  10ICE68  remains  below  10  km  and  quickly  reduces  to  values  smaller  than 
6  km  after  12  h.  Despite  this  small  tilt  magnitude  a  southerly  tilt  direction  is  preferred  until  late  in 
the  experiment.  Although  the  evolution  of  the  tilt  magnitude  is  similar  in  loRMNgg,  the  tilt  direction 
in  10RMN68  exhibits  more  variability  after  30  h  when  the  tilt  magnitude  is  small  (not  shown). 

405  3.4  SBC  formation  and  downward  flux  of  low-/),,  air  into  the  inflow  layer 

A  metric  for  the  downward  flux  of  low-0e  air  into  the  inflow  layer  (DFX)  has  been  defined  in  RMN 10 
(their  Sec.  4.2.1)  as 

DFX  =  w_0' ,  (4) 

where  w  denotes  downward  vertical  motion  and  0'e  the  perturbation  from  the  azimuthal  mean  of 
410  9e.  DFX  is  evaluated  at  1.5  km,  the  approximate  top  of  the  axisymmetric  inflow  layer. 

The  time-averaged  structure  of  the  SBC,  represented  by  low-level  upward  motion,  and  the  dis¬ 
tribution  of  DFX  shortly  after  the  time  when  vertical  wind  shear  is  imposed  is  depicted  in  Fig.  5. 
In  this  early  stage  of  the  evolution,  the  general  structure  of  the  SBC  and  the  distribution  of  DFX  is 
consistent  in  all  experiments.  The  SBC  forms  in  the  downshear-right  quadrant  extending  radially 
415  outwards  to  150  km  -  200  km.  Small  to  moderate  DFX  values  occur  over  a  broad,  banded  region 
downwind  of  the  SBC,  from  the  downshear  quadrant  through  downshear-left,  and  into  the  upshear 
quadrant.  In  all  experiments,  azimuthal  vertical  cross-sections  (as  Fig.  18  in  RMN,  not  shown  here) 
indicate  that  the  associated  downdrafts  form  underneath  the  helical  updrafts  of  the  SBC. 

It  is  of  interest  to  note  that  the  SBCs  in  RMN54,  ICE(;8,  and  njICEg*,  respectively,  extend  further 
420  into  the  upshear  quadrant  than  their  counterparts  in  RMN68,  CBLASTgg,  and  CBLAST54.  In  the 
former  experiments,  the  tilt  direction  at  early  times  exhibits  a  greater  downshear  (westerly)  com¬ 
ponent  than  in  the  latter  experiments  (Fig.  4).  The  low-level  vorticity  anomaly  associated  with  the 
tilted  vortex  can  thus  be  expected  to  extend  into  the  upshear  quadrant.  RMN  10  have  shown  evidence 
that  frictional  convergence  due  to  the  low-level  vorticity  anomaly  associated  with  the  tilted  vortex 
425  contributes  to  the  forcing  of  the  SBC.  The  correlation  of  the  SBC  formation  location  and  the  tilt 
direction  in  our  experiments  presented  here  provides  further  support  for  RMNlO’s  interpretation. 

A  distinct  maximum  of  DFX  is  found  in  the  upshear  quadrant  close  to  the  eyewall.  This  maximum 
is  most  pronounced  in  the  ICE  experiments  (Fig.  5(e)  and  5(f)).  In  RMN10,  this  DFX  maximum 
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was  not  given  special  attention.  A  pronounced  maximum  in  the  upshear  quadrant  occured  in  the 
430  20mps  case  of  RMN10  (their  Fig.  10(d)).  RMN10  further  noted  a  localized  region  of  downdrafts  in 
the  upshear  region  close  to  the  eyewall  after  the  TCs  in  the  lOmps  and  15mps  cases  have  re-aligned 
(their  Fig.  1 1(a)  and  associated  discussion). 

The  DFX  maximum  is  located  at  the  downwind  end  of  the  banded  region  of  positive  DFX  values. 
A  pronounced  0,  depression  is  found  at  the  top  of  the  inflow  layer  in  this  region  (not  shown).  Large, 
435  negative  values  of  0'e  therefore  contribute  to  high  values  of  DFX  (Eq.  4).  Downward  motion  into 
the  inflow  layer,  however,  is  particularly  strong  in  this  region  also  (not  shown).  A  preliminary  ex¬ 
amination  of  the  distributions  of  vertical  motion,  hydrometeors,  and  9e  suggests  that  the  downdrafts 
close  to  the  eyewall  may  form  by  evaporation  of  precipitation  falling  out  of  the  tilted  eyewall.  To  us, 
however,  it  is  not  clear  whether  the  more  pronounced  localized  maximum  in  the  ICE  experiments  is 
440  solely  due  to  the  inclusion  of  ice  microphysics  or  to  what  extent  the  differences  in  the  radial  structure 
of  the  TC  vortex  contribute.  Clearly,  more  research  is  needed  to  clarify  the  formation  mechanism  of 
these  downdrafts  and  the  relative  importance  of  the  confined  DFX  maximum  and  the  broad  region 
of  moderate  DFX  for  TC  intensity  evolution. 

3.5  Vortex-tilt  induced  dynamic  instability? 

445  A  non-linear,  Kelvin-Helmholtz  type  instability  in  the  inner  core  of  a  tilted,  dry  TC-like  vortex  has 
been  found  in  a  numerical  experiment  by  Reasor  et  al.  (2004,  their  Sec.  5a).  In  their  experiment, 
this  instability  contributed  to  enhanced  localized,  three-dimensional  mixing.  Such  mixing  is  one 
pathway  to  dissipate  kinetic  energy  in  the  TC’s  inner  core.  The  role  of  this  mixing  for  the  intensity 
evolution  of  TCs  in  vertical  wind  shear  is  hitherto  unexplored. 

450  The  necessary  condition  for  such  an  instability  is  that  the  gradient  Richardson  number,  1Z,  drops 
below  unity  (Abarbanel  et  al.,  1986).  Here,  7Z  is  calculated  as 

7 Z  =  N2/((du/dz)2  +  (dv/dz)2),  (5) 

where  N  is  the  (dry)  vertical  static  stability,  and  u  and  v  are  the  radial  and  tangential  wind  compo¬ 
nents10,  respectively.  The  tilt  of  the  vortex  may  lead  to  a  high  local  vertical  shear  of  the  swirling 
455  winds  in  the  inner  core.  Thus,  in  such  regions,  the  gradient  Richardson  number  may  reach  small 
values. 

In  all  of  our  (moist)  TC  experiments,  the  potential  for  a  Kelvin-Helmholtz  type  instability  and 
the  associated  three-dimensional,  turbulent  mixing  in  the  inner  core  is  an  element  of  vertical-shear 
interaction.  Underneath  the  eyewall  and  the  SBC,  TZ  <  1  is  an  ubiquitous  feature  within  the  layer  of 
460  strong,  frictionally-driven  inflow  (not  shown).  Above  this  inflow  layer,  localized  regions  of  TZ,  <  1 
occur,  at  least  temporarily,  in  all  experiments  presented  in  this  study  within  the  first  6  h  after  shear 

'^Normally,  the  definition  of  R  involves  the  zonal  and  meridional  wind  components.  Our  calculation  of  R  is  equivalent  to 
this  definition. 
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is  imposed.  These  regions  are  located  within  the  eyewall  in  the  downshear-left  to  upshear  quadrant 
(exemplified  in  Fig.  6  for  RMNfl8  and  ICEf,«).  This  feature  is  most  pronounced  in  RiVIN'(l,s  and  least 
pronounced  in  RMN54,  as  might  be  expected  for  the  weaker  vortex  in  RMN54.  The  regions  of  1Z  <  1 
465  are  frequently  located  radially  inwards  (within  20  km  -  30  km)  of  regions  of  high  DFX  (exemplified 
in  Fig.  6).  Due  to  this  proximity,  it  seems  justified  to  hypothesize  that  enhanced  mixing  associated 
with  the  tilt-induced  instability  locally  modifies  the  downward  flux  of  low-d,.  air  into  the  inflow 
layer.  A  quantitative  diagnostic  of  this  mixing  process  is  deferred  to  future  work. 

3.6  Inflow  layer  6e  depression 

470  3.6.1  0e  depression  outside  of  the  eyewall 

A  snapshot  of  the  height-averaged  inflow  layer  0e  distribution  after  5  h  into  the  shear  experiments 
is  presented  in  Fig.  7.  As  can  be  expected  from  the  DFX  patterns  discussed  above,  the  inflow  layer 
0e  distribution  is  similar  in  all  experiments  at  this  time:  Low -0e  air  wraps  cyclonically  from  the 
downshear  quadrant  into  the  upshear  quadrant.  Also  evident  in  all  cases  is  the  pronounced  wave- 
475  number  1  signature  of  the  SBC  in  the  vertical  motion  field  at  low  levels  above  the  inflow  layer. 

To  facilitate  the  intercomparison  of  the  experiments  we  consider  now  the  deviation  from  the  ax- 
isymmetric  mean  0,  distribution  at  the  time  when  shear  is  imposed  in  the  respective  experiment 
(Fig.  8).  This  figure  depicts  time-averaged  values  from  4  h  -  9  h  after  shear  is  imposed.  In  all  cases 
the  inflow  layer  0e  values  outside  of  the  inner-core  updrafts  are  significantly  depressed.  A  depres- 
480  sion  of  up  to  15  K  can  be  found  in  RMN68,  RMN54,  and  ICE68.  Minimum  (),.  values  are  generally 
found  to  the  southeast  of  the  centre  (upshear  to  downshear-left).  This  location  is  consistent  with  the 
DFX  maximum  in  this  region  and  the  continuous  depression  of  6e  values  along  the  spiralling  inflow 
leg,  as  hypothesised  by  RMN10  (their  Fig.  11a).  The  inflow  layer  0e  depression  extends  into  the 
downshear-right  quadrant  also,  outwards  to  a  radius  of  100  km  -  150  km.  In  the  absence  of  down- 
485  ward  flux  of  low  0e  in  this  region  (cf.  Fig.  5)  the  inflow  layer  0e  values  gradually  recover  while  air 
parcels  circle  from  the  upshear  quadrant  to  the  upshear-right  quadrant. 

The  inflow  layer  6e  depression  outside  of  the  eyewall  in  the  CBLAST  experiments  is  less  pro¬ 
nounced  than  in  the  comparable  RMN  experiments.  While  the  minima  in  RMN  exceed  -13  K,  the 
minima  in  CBLAST  have  values  of  approx.  -10  K.  This  less  pronounced  0e  depression  is  consis- 
490  tent  with  smaller  average  DFX  in  CBLAST(j8  and  CBLAST.54  as  compared  to  RMN@8  and  RMN54, 
respectively  (Fig.  5).  The  smaller  DFX  values  in  CBLAST  are  associated  with  a  weaker  general 
downdraft  activity  (see  section  3.7  and  Fig.  10  for  supporting  evidence). 

In  the  ICE  experiments,  the  spatial  pattern  of  the  0e  depression  is  radially  more  confined  than  in 
the  warm-rain  experiments.  This  is  consistent  with  the  DFX  pattern.  As  discussed  at  the  end  of 
495  Sec.  3.4,  it  is  not  evident  to  us  to  what  extent  this  greater  radial  confinement  can  be  attributed  to 
differences  in  the  cloud  microphysics  or  to  differences  in  the  radial  structure  of  the  TC  vortices. 
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3.6.2  9e  depression  underneath  the  inner-core  updrafts 


An  axisymmetric  view  of  the  inflow  layer  and  eyewall  9e  distribution  is  presented  in  RMN10.  In  this 
axisymmetric  view  a  tongue  of  low-//  air  effectively  intrudes  into  the  inflow  layer  from  above  and 
500  radially  inwards  towards  the  eyewall.  The  9e  values  of  air  rising  in  the  eyewall  of  the  sheared  TCs  are 
several  degrees  smaller  than  in  the  experiment  without  vertical  shear  (Fig.  9  in  RMN10).  The  same 
general  characteristics  of  the  axisymmetric  inflow  layer  and  adjacent  eyewall  structure  are  found  in 
the  additional  experiments  considered  in  this  study  (not  shown).  When  viewed  from  the  simplified 
perspective  of  an  idealised  Carnot-cycle  a  weakening  of  the  sheared  TCs  can  thus  be  expected.  As 
505  noted  in  RMN10,  however,  it  is  an  open  question  whether  an  axisymmetric  viewpoint  is  sufficient  in 
the  presence  of  pronounced  wave-number  1  asymmetries  in  the  inflow  layer  9e  distribution  and  the 
eyewall  updrafts  of  a  TC  in  vertical  wind  shear. 

In  the  present  study  we  complement  the  axisymmetric  view  by  the  examination  of  the  horizontal 
distribution  of  inflow  layer  9e  and  low-level  eyewall  updrafts.  For  the  following  discussion  the 
510  eyewall  updrafts  are  defined  as  the  contiguous  region  within  the  1  ms-1  vertical  velocity  isotach 
near  the  centre  in  Fig.  8.  We  consider  the  9e  values  within  this  region  to  be  a  good  approximation 
for  the  9e  values  of  air  rising  out  of  the  boundary/  sub-cloud  layer  in  the  eyewall  updrafts.  Upon 
examining  the  results,  it  is  obvious  that  in  all  cases  the  inner-core  updrafts  exhibit  a  pronounced 
wave-number  1  asymmetry  and  that  the  eyewall  is  not  closed. 

515  By  visual  inspection  we  estimate  the  average  9e  depression  underneath  the  eyewall  (averaged  from 
4  h  -  9  h)  to  be  0  K  for  CBLAST54,  2  K  -  3  K  for  RMN68  and  CBLAST68,  3  K  -  5  K  for  RMN54 
and  10ICE68,  and  6  K-  7  K  for  ICEg8.  These  values  correlate  well  with  the  intensity  evolution  of  the 
respective  TCs.  The  TC  in  CBLASTg4  does  not  weaken  until  7  h  after  shear  is  imposed  (Fig.  2(c)). 
The  TCs  in  RMN(l*  and  CBLASTg8  show  a  very  similar  weakening  in  the  first  12  h  while  R1VIN51 
520  and  yTCEgg  exhibit  slightly  faster  weakening  rates.  The  pronounced  and  rapid  weakening  of  the 
TC  in  ICEg§  is  associated  with  the  most  pronounced  depression  of  9e  underneath  the  eyewall.  This 
correlation  supports  RMNlO’s  result  that  the  intensity  evolution  of  resilient  TCs  in  vertical  shear  is 
governed  to  first  order  by  the  systematic  depression  of  9e  in  the  inflow  layer. 

The  inflow  layer  9e  depression  outside  of  the  eyewall  in  CBLASTg8  is  less  pronounced  than  in 
525  RMNgg.  The  9e  depression  underneath  the  eyewall,  however,  is  very  similar  in  both  experiments. 
Consistently,  the  intensity  decrease  is  very  similar  in  CBLASTgg  and  RMN68  in  the  first  12  h. 
It  can  therefore  be  inferred  that  the  replenishment  of  inflow  air  spiralling  towards  the  eyewall  is 
more  complete  in  RMN  than  in  CBLAST.  As  discussed  in  the  Introduction,  the  smaller  value  of 
Cd  in  CBLAST  may  lead  to  weaker  frictional  inflow  as  compared  to  RMN  and  thus  to  a  longer 
530  replenishment  time.  Here,  this  effect  of  a  smaller  Cjj  value  is  apparently  overruled  by  the  reduced 
extraction  of  heat  from  the  ocean  surface  due  to  a  smaller  value  of  Ck ■  It  is  difficult  to  examine  this 
aspect  in  RMN,54  and  CBLAST,54  because  the  9e  depression  underneath  the  eyewall  is  different  at 
this  early  time  in  the  experiment. 
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The  replenishment  time  of  air  with  depressed  9e  values  may  depend  considerably  on  the  azimuthal 
535  location.  The  radial  velocity  averaged  over  the  lowest  1  km  exhibits  a  pronounced  wave-number  1 
asymmetry  in  all  experiments  (Fig.  8).  In  the  downshear  semicircle  the  maximum  height-averaged 
radial  inflow  is  approx.  20  ms-1  while  in  the  upshear  semicircle  the  height-averaged  radial  velocity 
is  very  small  or  is  even  directed  outwards  (e.g.  in  RMN54  (Fig.  8(c))).  Figure  8  indicates  that 
the  bulk  of  the  depressed  inflow  layer  9e  spirals  around  the  centre  before  entering  the  inner  core 
540  updrafts  in  the  downshear-right  quadrant  (e.g.  Figs.  8(a)  and  8(f)).  In  RMN54  and  ICEgg,  in  contrast, 
it  is  indicated  that  low-(9e  air  from  the  downshear  quadrant  enters  the  eyewall  updrafts  on  a  much 
more  direct  pathway  in  the  downshear  to  downshear-left  quadrant  also  (Figs.  8(c)  and  8(e)).  The 
replenishment  of  inflow  layer  9e  along  this  more  direct  pathway  can  be  expected  to  be  less  complete 
than  along  the  longer  inflow  spiral.  In  the  presence  of  a  pronounced  wave-number  1  asymmetry  of 
545  the  inflow,  it  stands  to  reason  that  the  (azimuthal)  location  of  the  flushing  of  the  inflow  layer  with 
low-fi,  air  may  be  of  importance  for  the  intensity  evolution  also,  and  not  only  the  magnitude  of  the 
9e  depression  in  the  azimuthal  mean. 

3.7  Vertical  motion  characteristics 

The  formation  of  the  SBC  and  the  associated  downdrafts  that  flush  low-6*  air  into  the  inflow  layer 
550  are  two  central  aspects  of  RMNlO’s  framework.  We  here  examine  the  general  characteristics  of 
vertical  motion  in  our  suite  of  experiments  using  so-called  contoured  frequency  by  altitude  diagrams 
(CFADs  Yuter  and  Houze,  1995).  These  diagrams  depict  the  frequency  of  the  occurrence  of  vertical 
motion  of  a  given  magnitude,  here  binned  in  0.5  ms-1  intervals,  within  a  given  region  and  at  a  given 
altitude.  It  is  very  difficult  to  focus  specifically  on  the  SBC  and  on  regions  where  shear-induced 
555  downdrafts  form.  We  therefore  consider  all  vertical  motion  within  a  radius  of  160  km. 

A  comparison  of  the  RMN  and  ICE  experiments  is  presented  in  Fig.  9.  Without  vertical  wind  shear 
(Figs.  9(a)  and  9(b))  strong  up-  and  downdrafts  occur  more  frequently  in  the  ICE  experiment.  Under 
the  influence  of  vertical  shear  (Figs.  9(c)  and  9(d)),  the  frequency  distribution  broadens  considerably 
below  5  km  -  6  km  as  compared  to  the  respective  no-shear  experiment.  For  both,  the  warm  rain 
560  and  ice  microphysics  configurations,  vertical  wind  shear  increases  the  frequency  of  the  strongest  up- 
and  downdrafts  at  midlevels  and  below.  At  upper  levels  the  increase  of  downdraft  frequency  is  less 
pronounced  and  the  frequency  of  strong  updrafts  tends  to  decrease,  except  for  the  strongest  updrafts 
in  ICEgg.  It  is  evident  also  that  strong  downdrafts  below  5  km  occur  considerably  more  frequently  in 
ICEgg  than  in  RMNgg  (cf.  Figs.  9(c)  and  9(d)).  The  more  pronounced  downdraft  activity  in  the  ICE 
565  experiments  is  consistent  with  the  tenet  that  the  cooling  associated  with  the  melting  of  ice  particles 
may  initiate  or  strengthen  existing  downdrafts  (Srivastava,  1987). 

In  CBLAST,  the  imposed  vertical  wind  shear  leads  to  a  broadening  of  the  frequency  distribution 
also  (not  shown).  The  differences  between  CBLAST  and  RMN  experiments  are  less  pronounced 
than  between  ICE  and  RMN.  The  differences  are  depicted  in  Fig.  10  as  the  ratio  of  the  frequen- 
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570  cies  in  RMN68  and  CBLAST68  (Fig.  10(a))  and  RMN54  and  CBLAST54  (Fig.  10(b)).  The  results 
are  very  similar  for  RMN  and  CBLAST  without  imposed  vertical  shear  (not  shown).  Apparently, 
the  frequency  of  the  strongest  downdrafts  is  more  than  275  %  greater  in  RMN  than  in  CBLAST. 
Strong  updrafts  tend  to  be  more  frequent  in  RMN  also,  in  particular  so  in  RMN54  as  compared  to 
CBLAST54. 

575  Differences  in  the  thermodynamic  properties  of  the  ambient  air  in  which  downdrafts  occur  may 
contribute  to  the  observed  differences  in  downdraft  activity.  The  strength  of  a  downdraft  depends, 
inter  alia,  on  the  buoyancy  difference  between  the  descending  air  parcel  and  the  ambient  air  along  the 
parcel’s  downward  path.  A  vertically  integrated  measure  of  this  buoyancy  difference  is  the  downdraft 
convective  available  potential  energy  (DCAPE,  Emanuel,  1994,  his  Sec.  6.3).  The  derivation  of 
580  DCAPE  assumes  the  following  two-stage  process.  First,  evaporation  of  precipitation  cools  an  air 
parcel  isobarically  through  a  wet-bulb  process,  i.e.  until  saturation  is  reached.  Then,  the  parcel 
descends  pseudo-adiabatically  with  just  enough  evaporation  to  keep  the  parcel  saturated. 

We  have  calculated  DCAPE  using  the  azimuthal  mean  soundings  at  the  time  when  shear  is  im¬ 
posed.  Because  it  is  not  entirely  clear  which  radius  is  most  appropriate  to  examine  the  downdrafts 
585  underneath  the  SBC  we  considered  a  number  of  different  radii  and  radial  averages  between  60  km 
and  150  km  (not  shown).  For  all  of  these  regions,  a  clear  relationship  between  the  DCAPE  values 
and  the  downdraft  activity,  respectively  the  downward  flux  of  low-0e  air  into  the  inflow  layer,  is  not 
found.  In  particular,  DCAPE  values  below  500  hPa  in  CBLAST54  are  higher  than  in  RMN54,  in 
contrast  to  the  more  pronounced  downdraft  activity  in  RMN54. 

590  Apparently,  consideration  of  solely  thermodynamic  properties  of  the  air  in  which  the  downdrafts 
occur  is  insufficient  to  explain  the  observed  differences  in  downdraft  activity.  We  propose  the  follow¬ 
ing  alternative  explanation  that  focuses  on  the  dynamic  coupling  of  the  convection  with  the  frictional 
boundary  layer.  First,  we  note  that  the  weaker  downdrafts  in  CBLAST68  and  CBLAST54  are  asso¬ 
ciated  with  weaker  updraft  activity.  This  may  be  expected  because  a  prominent  downdraft  formation 
595  mechanism  is  the  evaporation  of  precipitation  that  is  formed  in  the  updrafts  and  falls  into  unsaturated 
air.  Both  the  smaller  value  of  Co  and  Ck  in  CBLAST  may  contribute  to  weaker  updrafts.  Ascent 
out  of  the  boundary  layer  forced  by  frictional  convergence,  either  underneath  the  eyewall  or  under¬ 
neath  the  vorticity  anomaly  associated  with  the  tilt  of  the  outer  vortex  (RMN  10,  their  Fig.  17),  can  be 
expected  to  decrease  for  decreasing  Cd-  F°r  example,  in  simple  Ekrnan  pumping  theory  the  vertical 
600  motion  at  the  top  of  the  boundary  layer  is  proportional  to  \J Co  (e.g.  Holton,  2004,  his  Eq.  5.38). 
The  smaller  value  of  Ck  in  CBLAST  decreases  the  rate  by  which  heat  is  extracted  from  the  ocean 
surface.  These  reduced  surface  fluxes  destabilise  the  boundary  layer  in  CBLAST  to  a  lesser  extent 
than  in  RMN.  Less  vigorous  updrafts  can  be  expected  to  occur  in  such  a  more  stable  stratification. 
Our  foregoing  argument  is  supported  by  the  results  of  the  CBLAST  and  RMN  no-shear  experiments 
605  (not  shown).  Without  vertical  shear,  the  activity  of  transient  rain  bands  and  the  associated  localized 
downward  flux  of  low  9e  into  the  inflow  layer  in  CBLAST  are  also  less  pronounced  than  in  the 
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no-shear  RMN  experiment. 

In  all  experiments,  vertical  wind  shear  increases  the  frequency  of  moderate  to  strong  downdrafts. 
It  is  clear  that  CBLAST  has  the  lowest  frequency  of  strong  downdrafts  and  that  ICE  has  the  highest 
610  frequency.  The  downdraft  activity  as  seen  in  the  CFADs  is  consistent  with  the  magnitude  of  the 
inflow  layer  9e  depression  outside  of  the  eyewall  (cf.  Fig. 8).  Therefore,  the  examination  of  the 
CFADs  provides  additional  support  for  the  hypothesis  that  the  distinct  depression  of  the  inflow  layer 
6e  values  are  caused  by  downdrafts  induced  by  the  imposed  vertical  wind  shear. 

4  A  closer  look  at  the  ICE6s  experiment 

615  The  weakening  of  the  TC  in  ICEgg  is  much  more  pronounced  than  in  the  comparable  warm-rain 
experiments  RMN(iS  and  CBLASTfj8.  In  the  previous  sections  it  has  been  shown  that  the  general 
structural  changes  in  the  first  20  h  are  nevertheless  very  similar  in  all  experiments.  We  now  take  a 
closer  look  at  ICE(i8  to  further  confirm  that  the  intensity  evolution  in  this  case  is  consistent  with  the 
shear-induced  thermodynamic  modification  of  the  inflow  layer.  First,  the  inflow  layer  9e  distribution 
620  is  examined  during  periods  of  distinct  intensity  changes.  Then,  we  will  provide  evidence  based 
on  axisymmetric  spindown  theory  that  the  rapid  weakening  in  the  first  20  h  is  consistent  with  the 
spindown  of  the  TC  vortex  due  to  a  significant  decrease  of  the  inner  core  convective  mass  flux. 

4.1  Intensity  evolution  and  inflow  layer  9e 

It  has  been  shown  in  section  3.6  (Fig.  7(e))  that  the  rapid  weakening  is  associated  with  a  significant 
625  depression  of  inflow  layer  9e  soon  after  shear  is  imposed.  Here  we  emphasize  that  a  significant  part 
of  the  depressed  0e-air  apparently  intrudes  from  the  downshear  quadrant  into  the  inner-core  updrafts 
on  a  direct  pathway  (Fig.  11(a),  at  7  h).  Over  half  of  the  area  underneath  the  inner-core  updrafts  at 
this  time  has  9e  values  lower  than  355  K,  as  compared  to  approx.  363  K  before  shear  is  imposed 
(not  shown).  Along  with  this  significant  frustration  of  the  9,  values  in  the  inner-core  region  the 
630  vertical  mass  flux  at  midlevels  (5  km)  decreases  in  the  first  10  h  by  60  %  from  50  x  106kg  s-1  to 
approx.  20  x  106kg  s-1  (not  shown). 

We  now  consider  the  two  distinct,  short-lived  local  intensity  maxima  at  22  h  and  44  h,  prominent 
features  of  the  intensity  evolution  in  ICE(is  (Fig.  2(b),  the  maxima  are  found  at  64  h  and  86  h 
because  this  figure  includes  the  42  h  spinup  period).  Both  periods  of  temporary  re-intensification  are 
635  preceded  by  relatively  high  inflow  layer  9e  values  in  the  downshear  to  downshear-right  quadrants 
while  the  very  pronounced  9e  depression  at  these  times  is  mostly  confined  to  the  upshear  semicircle 
(at  20  h  (Fig.  1 1(b))  and  at  40  h  (Fig.  11(d))).  Due  to  the  distinct  wave-number  1  inflow  asymmetry 
discussed  previously  (Fig.  8),  it  can  be  expected  that,  at  these  times,  the  eyewall  updrafts  are  fed  to 
a  large  extent  by  the  relatively  high-d,  air  in  the  downshear  semicircle.  Eyewall  0,  values  may  then 
640  temporarily  increase  and  the  TC  may  intensify.  The  temporary  intensification  ensues  until  the  very 
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low  9e  air  from  the  upshear  semicircle  wraps  around  the  centre  and  limits  considerably  the  supply 
of  relatively  high  9e  (exemplified  at  23  h  (Fig.  1 1(c))  and  at  45  h  (Fig.  1 1(e))). 

At  the  end  of  our  experiment  a  robust  re-intensification  is  indicated  in  the  intensity  time  series 
(Fig.  2(b)).  As  in  the  experiments  examined  in  RMN10,  this  final  re-intensification  is  associated 
645  with  an  almost  complete  recovery  of  the  inflow  layer  9e  (Fig.  1 1(f)). 

As  indicated  in  Figs.  1 1(b)  -  1 1(e),  low-0e  air  in  the  inflow  layer  may  wrap  around  the  centre  into 
the  downshear-right  quadrant.  This  low-((,.  air  tends  to  suppress  the  convection  of  the  SBC.  Further¬ 
more,  after  20  h  in  rCE(j«  the  vortex  tilt  is  no  longer  in  the  approximate  left-of-shear  equilibrium 
direction  but  starts  to  precess.  In  accordance  with  this  precession,  the  vorticity  anomaly  associated 
650  with  the  tilt  of  the  outer  vortex  then  starts  to  move  also  and  thus  no  longer  provides  persistent  forc¬ 
ing  for  the  SBC  in  the  downshear-right  quadrant.  It  is  interesting  to  note  that  a  narrow  band  of 
convection  preferably  occurs  in  the  downshear  quadrant  at  these  later  times  (Figs.  11(b)  -  11(e)). 
Frequently,  one  observes  significant  downward  flux  of  low  9e  associated  with  this  remnant  SBC  (not 
shown). 

655  4.2  A  timescale  for  vortex  spindown 
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Within  the  simplified  framework  of  an  idealised  Carnot  cycle,  and  based  on  the  idea  that  environ¬ 
mental  vertical  wind  shear  leads  to  the  frustration  of  such  a  thermodynamic  cycle,  Tang  and  Emanuel 
(2010)  have  derived  a  theory  for  the  storm’s  reduced  steady-state  intensity.  Such  a  theory,  however, 
does  not  provide  information  about  the  temporal  evolution  of  storm  intensity,  i.e.  how  rapidly  a  TC 
weakens.  A  well-known  concept  for  the  weakening  of  a  TC  is  the  axisymmetric,  frictionally-induced 
spindown  of  the  TC  vortex.  Frictional  convergence  forces  ascent  out  of  the  inflow  layer.  When  the 
vertical  convective  mass  flux  is  smaller  than  the  frictionally  forced  vertical  mass  flux,  some  of  the 
air  has  to  diverge  above  the  inflow  layer  due  to  mass  continuity  (Willoughby,  1979).  Under  the 
assumption  that  air  parcels  approximately  conserve  their  angular  momentum  above  the  frictional 
inflow  layer,  this  divergence  implies  a  decay  of  the  tangential  winds,  i.e.  the  TC  weakens. 

A  timescale  for  the  frictional  spindown  of  an  axisymmetric  vortex  without  convective  mass  flux 
has  been  derived  by  Elliassen  and  Lystad  (1977).  Their  model  assumes  a  neutral  stratification  above 
the  layer  of  strong  frictional  inflow  which  is  a  reasonable  assumption  in  the  inner  core  of  a  TC 
because  isosurfaces  of  9e  and  angular  momentum  are  approximately  congruent  in  this  region.  Mont¬ 
gomery  et  al.  (2001)  have  verified  by  numerical  experiments  that  the  theory  gives  reasonable  results 
for  the  spindown  of  TC-like  vortices,  even  beyond  the  formal  validity  of  the  theory.  The  half-time 
that f,  i.e.  the  time  that  is  required  to  reduce  the  initial  tangential  velocity  v  by  a  factor  of  2,  is  given 
by  Elliassen  and  Lystad’s  theory  as 


H-h 

thalf  ~yj2 


(6) 


X2Cdv 

Here,  H  is  the  fluid  depth,  h  the  height  of  the  boundary  layer,  and  %  a  reduction  factor  defined  by 
Elliassen  and  Lystad  as  the  ratio  of  the  tangential  winds  at  the  top  of  the  surface  layer  to  the  tangential 
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winds  at  the  top  of  the  Ekman  layer.  Estimating  H  —  h=  104  m,  \  =  0-8,  and  using  Cp  =  3.8  x  10~4 
from  Deacon’s  formula  with  v  =  68  ms-1  yields  a  half-time  thaif  =  16.8  h. 

The  weakening  of  the  TC  in  ICEgs  is  remarkable:  intensity  decreases  from  68  ms-1  at  42  h  to 
680  33  ms-1  at  68  h  (Fig.  2(b)),  corresponding  to  a  weakening  rate  of  35  ms'1  in  26  h;  or  a  half-time 

of  approx.  25  h.  This  value  is  consistent  with  a  partial  shutdown  of  the  convective  mass  flux.  As 
noted  above,  concurrently  with  the  frustration  of  the  inflow  layer  0,  and  the  rapid  weakening,  the 
vertical  mass  flux  at  5  km  height  decreases  by  approx.  60%.  Consequently,  it  can  be  argued  that  the 
significant  reduction  of  the  convective  vertical  mass  flux  by  the  severe  depression  of  the  inflow  layer 
685  9e  values  broadly  governs  the  rapid  weakening  in  ICE68. 

5  Conclusions 

In  the  light  of  a  new  framework  developed  in  a  recent  study  (RMN10),  the  current  study  further 
examines  the  evolution  of  TCs  exposed  to  environmental  vertical  wind  shear.  The  new  framework 
has  been  developed  based  on  idealized  numerical  experiments  for  which  we  have  deliberately  chosen 
690  a  simple  set  of  physical  parameterisations  to  minimize  the  complexity  of  the  problem.  In  the  current 
study,  we  survey  and  diagnose  five  additional  numerical  experiments  with  some  modifications  of  the 
experimental  setup  to  assess  the  robustness  of  our  previous  results. 

The  modifications  of  the  experimental  setup  focus  on  parameters  that  are  important  for  structure 
change  and  intensity  evolution  of  numerically  simulated  TCs  and  the  storms’  evolution  in  vertical 
695  wind  shear,  discussed  in  some  detail  in  the  Introduction  in  Sec.  1.2.  Specifically,  in  comparison 
to  RMN10,  we  here  consider  weaker,  but  still  mature  TCs;  employ  values  for  the  exchange  coeffi¬ 
cients  of  surface  heat  and  momentum  fluxes  that  are  in  better  agreement  with  recent  observational 
results;  and  complement  our  previous  warm-rain  experiments  with  experiments  that  employ  an  ice 
microphysics  scheme. 

700  The  ability  of  TC-like  vortices  to  withstand  the  differential  advection  by  vertical  wind  shear  and  to 

eventually  re-align  depends  on  the  radial  profile  of  the  swirling  winds  (discussed  in  detail  in  Sec. 3.1). 
Some  differences  in  this  radial  structure  arise  in  the  individual  experiments  before  interacting  with 
the  imposed  vertical  wind  shear.  The  most  notable  differences  are  found  in  the  experiment  with  ice 
microphysics  in  which  the  TC,  inter  alia,  features  a  radius  of  maximum  winds  that  is  approx.  25  % 
705  smaller  than  in  the  warm-rain  experiments.  Differences  in  the  radial  vortex  structure  complicate  the 
attribution  of  differences  in  the  TC  evolution  to,  e.g.,  the  inclusion  of  ice  microphysics  alone. 

The  differences  in  the  vortex  structure  notwithstanding,  in  all  experiments  the  TCs  exhibit  the 
same  general  structural  changes  when  exposed  to  vertical  wind  shear.  The  structural  changes  are 
consistent  with  RMN  10’s  new  framework  and  comprise  a  vortex  tilt  predominantly  to  the  downshear- 
710  left  for  20  h  -  24  h,  the  formation  of  a  pronounced  wave-number  1  convective  asymmetry  outside 
of  the  eyewall  (reminiscent  of  a  stationary  band  complex),  and  downward  flux  of  low-0e  air  into  the 
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inflow  layer  in  the  downshear-left  semicircle.  A  significant  depression  of  the  inflow  layer  9e  values 
(0(10  K ))  occurs  in  all  experiments. 

Furthermore,  for  all  experiments  the  potential  for  a  non-linear,  Kelvin-Helmholtz  type  instabil- 
715  ity  in  the  inner  core  of  the  tilted  TC-vortex  is  diagnosed.  Such  an  instability  would  contribute  to 
enhanced  dissipation  of  kinetic  energy  and  may  modify  the  flux  of  low-fi,  air  into  the  inflow  layer. 
A  quantitative  analysis  of  the  mixing  process  associated  with  the  potential  instability,  however,  is 
deferred  to  future  work. 

Two  additional  diagnostics  are  considered  in  the  current  study  that  have  not  been  considered  in 
720  RMN10,  namely  the  0,  values  underneath  the  asymmetric  inner-core  updrafts,  and  so-called  con¬ 
toured  frequency  by  altitude  diagrams  to  examine  the  general  characteristics  of  vertical  motion.  The 
magnitude  of  the  9e  depression  underneath  the  asymmetric  inner-core  updrafts  early  in  the  respec¬ 
tive  shear  experiment  correlates  well  with  the  ensuing  weakening  of  the  TC,  supporting  the  idea 
that  it  is  the  depression  of  the  inflow  layer  9e  that  governs  the  intensity  evolution  to  first  order.  For 
725  all  experiments,  vertical  wind  shear  increases  the  frequency  of  strong  downdrafts  as  compared  to 
the  no-shear  experiments.  The  magnitude  of  the  9e  depression  outside  of  the  eyewall  is  consistent 
with  the  frequency  of  strong  downdrafts,  supporting  the  idea  that  it  is  the  shear-induced  downdrafts 
that  flush  low-0e  air  into  the  inflow  layer.  Consequently,  both  additional  diagnostics  corroborate 
RMNlO’s  framework. 

730  In  one  of  our  experiments  employing  ice  microphysics,  a  very  rapid  weakening  of  the  TC  oc¬ 
curs.  This  intensity  evolution  is  distinct  from  the  intensity  evolution  in  our  warm-rain  experiments. 
Therefore,  a  question  may  naturally  arise  whether  ice  microphysics,  as  represented  in  the  RAMS 
model,  introduces  an  important  mechanism  for  intensity  change  that  is  not  considered  in  RMNlO’s 
framework.  As  noted  above,  the  structural  changes  caused  by  the  vertical  wind  shear  in  the  experi- 
735  ment  with  ice  microphysics  are  very  similar  as  in  the  warm-rain  experiments.  In  the  ice  experiment, 
the  general  downdraft  activity,  the  flux  of  low-(9e  air  into  the  inflow  layer,  and  the  0,  depression 
underneath  the  inner-core  updrafts  are  all  the  most  pronounced  in  our  suite  of  experiments.  The 
inner-core  9e  depression  that  accompanies  the  initial  rapid  weakening  is  more  than  twice  as  high 
than  in  the  comparable  experiments  employing  warm-rain  microphysics  (6  K  -  7  K  as  compared  to 
740  2  K  -  3  K).  Axisymmetric  spindown  theory  (Elliassen  and  Lystad,  1977;  Montgomery  et  al.,  2001) 

further  supports  the  hypothesis  that  the  rapid  initial  weakening  is  due  to  a  substantial  decrease  of 
the  convective  mass  flux  associated  with  the  significant  depression  of  the  inflow  layer  9e  by  shear- 
induced  downdrafts.  In  addition,  we  have  documented  a  close  relationship  between  the  evolution  of 
the  inflow  layer  9e  values  and  two  transient  intensity  maxima  that  occur  at  22  h  and  44  h  after  verti- 
745  cal  wind  shear  is  imposed.  The  final  re-intensification  of  the  TC  in  this  experiment,  as  in  RMN10, 
is  associated  with  a  near-complete  replenishment  of  the  inflow  layer  9e  values  also. 

Based  on  the  evidence  presented  here,  we  conclude  that  RMNlO’s  framework  for  the  intensity 
modification  of  TCs  in  vertical  wind  shear  is  robust  to  the  presented  changes  in  the  setup  of  the 
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numerical  experiment.  We  advocate  further  study  and  tests  of  this  framework  with  case  studies  and 
750  real  atmospheric  data. 
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Table  1.  Overview  of  the  experiments  that  are  considered  in  this  study.  The  experiments  RMN68  and  10RMN68 
are  described  in  detail  in  RMN10.  In  column  Ck.  ’Deacon’  denotes  that  Deacon’s  formula  (see  section  2.3) 
is  used  to  determine  the  exchange  coefficients.  The  last  column  gives  the  TC’s  intensification  rate  just  before 
vertical  shear  is  imposed. 


Experiment 

'initial' 

intensity  [ms-1] 

Ck 

Ck /Cd 

microphysics 

shear 

[ms_1/12  km] 

preceding  int. 

rate  [ms_1/h] 

RVINfis 

68 

Deacon 

1 

warm 

15 

1.8 

10RMN68 

68 

Deacon 

1 

warm 

10 

1.8 

RMN54 

54 

Deacon 

1 

warm 

15 

0.0 

CBLASTes 

68 

1.16  x  1CT3 

0.65 

warm 

15 

1.1 

CBLAST54 

54 

1.16  x  1CT3 

0.65 

warm 

15 

0.0 

ice68 

68 

Deacon 

1 

ice 

15 

2.5 

10ICE68 

68 

Deacon 

1 

ice 

10 

2.5 
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Fig.  1.  (a)  Radial  profiles  of  azimuthal  mean  vorticity  at  2  km  (thin)  and  azimuthal  mean  tangential  wind  at  the 
radius  of  maximum  winds  at  1  km  height  (thick)  at  the  time  when  vertical  wind  shear  is  imposed  in  RMN54 
(red),  CBLASTeg  (blue),  and  ICEgg  (grey),  respectively,  (b)  Radial  vorticity  profiles  as  in  (a)  but  zoomed  in  on 
the  “skirt”  region  between  50  km  and  100  km. 
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Fig.  2.  a)  Intensity  time  series  for  the  RMN  no-shear  experiment  (red),  RMN54  (light  red),  RMNgs  (dark  red), 
the  CBLAST  no-shear  experiment  (blue),  CBLAST54  (light  blue),  and  CBLASTgg  (dark  blue).  The  colour- 
coded  arrows  denote  the  time  when  vertical  wind  shear  is  imposed  in  the  respective  experiment,  b)  As  in 
a),  but  for  the  RMN  no-shear  experiment  (red),  RMNgg  (dark  red),  ioRMNgs  (light  red),  the  ICE  no-shear 
experiment  (black),  ICEgs  (dark  grey),  and  loICEgs  (light  grey),  c)  Time  series  of  the  difference  between  the 
intensity  in  RMN54  (light  red),  RMNgg  (dark  red),  CBLAST54  (light  blue),  and  CBLASTgs  (dark  blue)  and 
the  respective  no-shear  experiment.  Time  0  h  denotes  the  time  when  vertical  wind  shear  is  imposed  in  the 
respective  experiment,  d)  As  in  c),  but  for  RMNgs  (dark  red),  ioRMNgs  (red),  ICEgs  (dark  grey),  and  10ICE68 
(light  grey). 
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Fig.  3.  a)  Time  series  of  the  tilt  magnitude  for  RMN54  (light  red),  RMN68  (dark  red),  CBLAST54  (light  blue), 
and  CBLASTs8  (dark  blue).  Vortex  tilt  is  defined  as  the  vector  difference  between  the  vorticity  centroids  at 
10  km  and  1  km  height.  For  visual  clarity,  the  3  h  running  mean  of  the  time  series  is  shown,  b)  As  in  a)  but 
for  RMNg8  (dark  red),  mRMNss  (light  red),  ICE68  (dark  grey),  and  10ICE68  (light  grey).  Note  that  the  tilt 
magnitude  for  ICE68  goes  off  the  scale  between  33  h  and  38  h  reaching  a  maximum  value  of  195  km  (not 
shown). 
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Fig.  4.  Time  series  of  the  tilt  direction  for  RMNes  (a),  CBLAST68  (b),  RMN54  (c),  CBLASTs4  (d),  ICEgs  (e), 
and  10ICE68  (f).  Vortex  tilt  is  defined  as  the  vector  difference  between  the  vorticity  centroids  at  10  km  and  1  km 
height.  Dashed  line  denotes  the  3  h  running  mean. 
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Fig.  5.  Downward  flux  of  9'e ,  DFX  (colour,  in  0. 1  Rms  1 ,  see  text  for  definition),  at  the  top  of  the  axisymmetric 
inflow  layer,  defined  as  1.5  km,  and  low-level  upward  lotion  (thin  contour:  0.25  ms-1,  thick  contour:  1  ms-1, 
averaged  from  1 .25  km  -  2  km  height)  both  averaged  from  1  h  -  6  h  after  vertical  shear  is  imposed  for  a)  RMNgs, 
b)  CBLAST68,  c)  RMN54,  d)  CBLAST54,  e)  ICEgs,  and  f)  roICEes-  The  arrow  indicates  the  direction  of  the 
deep-layer  shear.  The  horizontal  scale  is  in  km 


Fig.  6.  Snapshots  of  the  gradient  Richardson  number  at  3  h  after  vertical  shear  is  imposed  (shaded  for  values 
below  1)  in  (a)  RMN68  and  (b)  ICEgs.  Upward  motion  is  depicted  in  black  contours  (thin  for  0.25  ms-1,  thick 
for  1  ms-1).  The  cyan  contour  denotes  regions  of  strong  downward  flux  of  low-0e  air  into  the  inflow  layer 
(DFX  at  1.5  km  height,  the  contour  depicts  the  4.5  Kms-1  isoline).  The  horizontal  scale  is  in  km. 
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Fig.  7.  As  Fig.  5,  but  for  a  snapshot  of  inflow  layeg$e  (colour,  averaged  over  lowest  1  km)  and  low-level 
updrafts  at  5  h  after  vertical  shear  is  imposed  for  a)  RMNgs,  b)  CBLASTgs,  c)  RMN54,  d)  CBLAST54,  e) 
ICE68,  and  f)  10ICE68.  The  depicted  time  is  representative  for  the  early  part  of  the  experiments. 


Fig.  8.  Depression  of  inflow  layer  9e  (averaged  over  the  lowest  1  km  and  from  4  h  -  9  h)  relative  to  the 
axisymmetric  mean  at  the  time  when  vertical  wind  sh§a)'  is  imposed  (colours).  Thin  contours  depict  the  radial 
flow  averaged  over  the  lowest  1  km  (every  4  ms-1,  starting  at  2  ms-1,  dashed  for  negative  values  (outflow)). 
Thick  contour  denotes  low-level  updrafts  (1  ms~  ,  averaged  from  1.25  km  -  2  km  height),  a)  RMNes,  b) 
CBLASTgg,  c)  RMN54,  d)  CBLAST54,  e)  ICE68,  and  f)  i0ICE68 


Fig.  9.  Contoured  frequency  by  altitude  diagrams  of  vertical  motion  within  160  km  radius  of  the  centre.  Vertical 
motion  is  binned  into  intervals  of  0.5  ms~  ,  centred  on  the  values  given  on  the  x-axis.  The  colour-coded  value 
denotes  the  logarithm  of  the  frequency.  With  each  contour  the  frequency  decreases  by  a  factor  of  1/e.  Before 
taking  the  logarithm  of  the  frequency  a  small  positive  value  (10_s  «  e-11'5)  has  been  added  to  ensure  that  the 
logarithm  is  well  defined  everywhere.  For  the  no-shear  cases  of  RMN  (a)  and  ICE  (b)  the  6  h  time  period  after 
vertical  wind  shear  is  imposed  in  the  respective  shear  runs  is  considered.  For  RMNCg  (c)  and  ICEgg  (d)  the 
times  4  h  -  9  h  are  considered. 
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Fig.  10.  As  Fig.  9,  but  for  the  ratio  of  vertical  motion  frequency  of  RMNos/CBLASTes  (a)  and 
RMN54/CBLAST54  (b). 
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Fig.  11.  As  Fig.  7,  but  for  ICE68  at  a)  7  h,  b)  20  h.  c)  2^h,  d)  40  h,  e)  45  h,  and  f)  64  h. 


